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Abstract 

This paper addresses the classical and discrete Euler-Lagrange equations for 
systems of n particles interacting quadratically in K''. By highlighting the role 
played by the center of mass of the particles, we solve the previous systems via 
the classical quadratic eigenvalue problem (QEP) and its discrete transcendental 
generalization. The roots of classical and discrete QEP being given, we state 
some conditional convergence results. Next, we focus especially on periodic 
and choreographic solutions and we provide some numerical experiments which 
confirm the convergence. 
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1. Introduction 

This paper seeks to continue the development of the theory for the discrete 
calculus of variations which was initiated by Cresson and al., see 0, 0]- It 
consists originally in replacing the derivative x(i) of the dynamic variable x(<) 
defined on [io,i/] with a 2N + 1 terms scale derivative 

N 

□ex(0= ^^it + je)x-j{t), te[to,tf], 7jeC (1) 

j=-N 

where Xj{t) denotes the characteristic function of [to,tf] n [to +je,tf +je], for 
some time delay s. 

We consider a lagrangian C of n particles in R'', where d denotes the "phys- 
ical" dimension. The principle of least action may be extended to the case of 
non-differentiable dynamic variables. For conservative systems, the equations 
of motion may be returned as the following two dynamic sets of equations 

Xj(i) = Fj(xi, . . . ,x„) and - n..^n^Xj{t) ^ Fj{xi, . . . (2) 
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where the functions Fj , Fj are built from the specific interaction between the 
particles. While the first system in ([2]) deals with ODE, the second one consists 
in a set of functional difference equations. We investigate for each system the 
existence of pseudo-periodic solutions of the shape 



where Ug and constitute a family oi K + 1 vectors of and {Xe)g G (C*)-^ 
is a sequence of K distinct complex numbers. 

The rest of this paper is organized as follows. Section 2 is devoted to the 
derivation of the classical and discrete Euler-Lagrange equations (respectively 
abbreviated as C.E.L. and D.E.L.) and highlights the role played by the center 
of mass ;^Xs(t) where Xs{t) = J^j^jW- Section 3, we present a method 
for solving the equations of motion for generic lagrangians for C.E.L. as well as 
D.E.L. . The first step of this method determines Xs(i) from some generalized 
(quadratic or transcendental) eigenvalue problem. The second step seeks Xj(i) 
from :x.g(t) by solving another eigenvalue problem. Section 4 is devoted to the 
convergence of the generalized eigenvalue problem linked to the D.E.L. as e 
tends to 0. Section 5 deals with the existence and the features of periodic and 
choreographic solutions. Finally, Section 6 presents some numerical experiments 
illustrating the phenomenon of convergence as e tends to 0, for some various 
operators O^. 

2. Equations of motion for symmetric quadratic lagrangians of n par- 
ticles systems in M"^ 

The principle of least discrete action has been developed in 0, Q to which 
we refer throughout the paper. We denote by Cpw the space of the functions 
X : [^Oj — ^ continuous on each interval [to+je, to + {j + l)e] n [to, for all 
j >0 and small enough, i.e. j < ^ ~ ° . If X = (xi, . . . x„) denotes a system of 
n functions in Cp„ , we may think of X as the set of dynamic variables describing 
the state of a system of n interacting particles in W^. 

We consider actions Aconti-^disc ■ — >■ M of the shape 



Ao„t(X)=/ C{X{t),X{t))dt, Ad^sc{X)^ C{X{t),n,X{t))dt. (4) 



From now on, we drop t from the formulas when it is clear enough. 

We introduce a general quadratic lagrangian of n particles in R'', compatible 
with discrete symmetries of the system. Let Ji,..., J5 € C°{[to,tf],R'^^'^), 
Ji, . . . , J4 be symmetric matrices and Jg, J7 G C°{[to,tf],R'^). For an isolated 
particle with position x and velocity y we may set 



K 




(3) 
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Next, two particles with positions x^-, x^, and velocities y^, are interacting 
for pairs in conformity with the following lagrangian 

'C2(xj,y-,-,Xfe,yfc) = Jayfc + '^^jj^yik- 
Therefore, the lagrangian of the whole system is 

n 

£ = ^£i(xj,Xj) +^/:2(xj,Xj,Xfc,Xfc). (5) 

3 = 1 j=ik 

n 

Theorem 2.1. Lei X = (xi, . . . , x„) in Cp^ and Xg — Xj . 

A necessary and sufficient condition for 1^ to he a critical point of Acont 'in 
is that X satisfies the dynamic system 

( Jl - 2 J3)Xj + (-2^3 + ^5 + Ji - J5)Xj + Ch- J2+ 2Ji)Xj = 

-2J3i, -2j3X, +2J4X, + (J7- jg), Vj. (6) 

A necessary and sufficient condition for to be a critical point of Adisc in Cp^ 
is that X satisfies the linear functional recurrence system of equations 

□ - 2J3)aXj) + U^.CJ^Xj) + JsDeXj + (J2 - 2J4)x, = 

-2n_,(j3nex,)-2j4X3-(n_,j6 + J7)), Vj. (7) 



Proof. We first recall the classical and discrete Euler-Lagrange equations, which 
are respectively given by 

d r)C ■ BC 

rf^a-(^'^W'^w) + Q^^{t,m,m) = (8) 

and 

□-e^^(i,x(t),neX(0) + ^(t,x(0,neX(i)) = 0, (9) 

for all j G {1, . . . , n}. 

The computation of gradients of C needs the following property : if a, b e M"*, 

Vx( "^ax + ■^xb) a + b. 
Let us prove Because of the symmetry of Ji, J2, J3, Ji, we get 

= JiXj + 2J3 ^Xfc + V5XJ + Je and 

dC 

- — = J2Xj + 2 J4 Xfc + J^kj + Jt- 
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Then, by setting Xg — Xj , we get 

-— = (Ji - 2 J3)xj + 2 J3XS + ^JsXj + Jg and 
axj 

a£ 

- — = ( J2 - 2 J4)xj + 2 J4XS + J5XJ + J7. 
axj 

The equation dH) gives for aU j : 

( ii - 2 j3)xj + ( Ji - 2 J3)xj + ^isXj + V5XJ + ie + 2 j3Xs + 2 J3XS = 
( J2 — 2 J4)xj + J5XJ + J7 + 2 J4XS 

which is equivalent to (jS]). 

The proof of ([7]) is quite similar since © gives for all j : 

□ _e(JineXj - 2J3neXj) + D^^l VgXj) + JsDeXj + (J2 - 2J4)Xj = 

-2n_,( J3n,x,) - 2 J4X, - Je + h)- 

which implies ([7]). □ 

We notice that the equations ([6]) and ([7|) are quite uncoupled since the 
coupling is realized only through the vector x^. We mention two simple conse- 
quences of the previous result. The first one arises from summing all equations 
([6]) or summing all equations in ([7]) over and the second one deals with time- 
independent lagrangians such that J5 is skew-symmetric. 

Corollary 2.1. IfK is a solution to (d^, then the sum x^ satisfies the dynamic 
system 

(Ji + 2(n - 1) J3)x; + (ii + 2{n - 1)^3 + V5 - J5)x, - 

(J2 + 2{n - 1) J4 ~ ^j5)x. -I- n( jg - J7) - 0. (10) 

Similarly, if K is a solution to then x^ satisfies the functional equation 

□ _e(JineX,) + 2{n - l)n_e(J3ax.) + D-el^^SXs) + 

JsDeXs + ( J2 + 2{n - l)J4)xs + Je + nJj = 0. (11) 



Corollary 2.2. Suppose that the functions Jk{t), k G {1,...,7}, are constant 
w.r.t. time and J5 is skew-symmetric. The systems of equations ilO\) and (0) 
simplify respectively into 

(Ji + 2{n - 1) J3)x; - 2 J5X, - ( J2 + 2{n ~ 1) J4)x, = nj-j (12) 

and 

( Ji - 2J3)xj - 2J5XJ - ( J2 - 2 J4)xj = -2 J3X, -f 2J4X, + J7. (13) 
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Similarly, the systems of functional recurrence equations Ul\) and ^ simplify 
respectively into 

(Jl + 2(n - 1) J3)n_eax, + J5(ax. - D-eXs) + (./2 + 2{n - 1) J4)x, = 

-nn_e(l) - nJr (14) 



and 



(Jl - 2J3)a^,a,Xj + J5{a,Xj - D-eXj) + (Ja - 2J4)xj = 

-2 JgD-eax. - 2 J4X. + D-e - Jt- (15) 



Remark 2.1. Let Js,J9,Jw e C"([io, t/], M""'^"'') denote the matrices con- 
structed by blocks as follows 



2J3 
V2J3 



2J3 
Jl 



2J3 



2J3\ 



2J3 
^1/ 



Jo = 



2J4 
V2J4 



2J4 
J2 



2J4\ 



■ • . 2 J4 

2J4 J2 / 



and Jio = diag{J5, . . . , J^). If X e Cp^ is a critical point of Acont, i-e. 
satisfies ([6|) and if Y e C^^ vanishes at and tj, then we have 



A-contO^ + Y) — Aco7it0^) — 2 



( Y JsY + JgY + ^ JioY)di. 



to 



As a consequence, if the integrand is a positive definite quadratic form w.r.t. 
(Y, Y), then the equations © are necessary and sufficient conditions for a strict 
minimum of the action Acont to occur. Especially, if J5 = and the matrices 
Js and Jg are definite positive, such an optimum occurs. 



Remark 2.2. Let A € R''-^'^ a nonsingular matrix and b e 
consider the transformation of the whole system 

x,(i) = Ax,(t)+b. 



be given. Let us 



Then this transformation is covariant for quadratic lagrangians in the sense that 
£(xi, . . . , x„) is of the shape JS]) iff >C(xi, . . . , x„) is of the shape JS]). Moreover 
the properties of symmetry for Ji, . . . , J5 are equivalent to those for Ji, . . . , J5. 
At last, the equations of motion ([2]) are covariant altogether as can be shown 
from the formula for affine forces 



F,{±) = AF,{A-\^^ - b), . . . , A-i(x„ 
where X = (xi, . . . ,x„). 



b)). 
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Remark 2.3. If J5 = 0, the system has a lagrangian of the shape T(X) — f/(X) 
and consequently, it is conservative, i.e. the energy 

" /I 1 \ 

"^J-^i^J ~ 9 ^^i-^axj + ^Jeij - ^JtXj + ^ (^Xj Jgifc - ^Xj J4Xfc) 

(16) 

is a constant of motion. 



3. Solutions of equations of motion in the general case 

3.1. Preliminaries on Quadratic Eigenvalue Problems 

We provide in this section the solutions to problems presented in Corollaries 
12.11 and 12.21 From now on we suppose that the vectors and matrices Jfc , k = 
1, . . . , 7, are time- independent and that J5 is skew-symmetric. By general case, 
we mean that the set of coefficients ( Jfc)fc=i,...,7, satisfying the conditions (fT9]) 
and ^ below is everywhere dense in (R'^'"^)^ x (R'')^. 

According to define the Quadratic Eigenvalue Problem associated 

to [A, B, C) as the search of the complex roots of the discriminantal equation 

det(AA2 +BX + C)=0 (17) 

where the l.h.s. is a polynomial of A of degree 2d, together with the description 
of the various kernels ker(AA^ -|- BX + C). The reader is referred to 9] for a 
survey of theory applications and algorithms of the QEP. 

It is a classical old fact 0, i, i] that, if all the roots A^, £ G {1, . . . , 2d}, of 
A\^-\-BX-\-C are distinct and C is invertible, the general solution to the dynamic 
system + Sx + Cx ~ ko has the shape ([3]), with K = 2d, Uq C^^ko and 
Ml e ker(A^A -I- XgB -t- C). When the number of roots is less than 2d, slight 
more complicated expressions for the solutions may be found in [^, Q 
Because of equations (fTO)) to p3|) and the previous discussion, we may provide, 
under specific assumptions, the shape of the solutions to C.E.L. and D.E.L.. 

3.2. The case of C.E.L. 

We introduce the matrix-valued function 

V,{X) (Ji + 2(j/ - 1)J3)X^ - 2J5A - (J2 + 2(1/ - 1) J4) (18) 
and the following subsets of C 

:= {X e C/det{V^{X)) = 0}, Vi^ G K. 
Proposition 3.1. We assume that 

\Qn\ = IQol = 2d, Q„nQo-0, 

det(-p„(0)) ^ and det(7'o(0)) ^ (19) 

where n denotes the number of interacting particles. Then all the solutions Xg 
and Xj to and hl!^) respectively are of the shape ^ with K — 2d. 
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Proof. Since det(J2 + 2(n — 1) J4) 7^ 0, we may define 

Xs,o = -n{J2 + 2{n - l)J^y^Jj 

and we have obviously ^ Q„. The first condition |Q„| = 2d guarantees that 
the solution Xs to is of the shape ^ with K = 2d, i.e. 

x,(t) =x,,o+ e"*x,,„, (20) 

for some convenient vectors Xs,a S ker (7^„(a)). Since Q„ H Qo = 0, the matrix 
Vo{ci) is invertible for each a G Qn so we may define for j G {1, . . . , n} 

Let us set Xj^o = — (J2 — 2 J4)^^(2 J4XS.0 + "Z?)- Straightforward computations 
show that Xj — -^^8,0 and Xj = 'k^s,a, Va G Q„. Therefore, we have proved 
that Xj +X)«eQ — n^s(^) ^ particular solution to (IT^ . Hence, the 

general solution to is given by the formula ^ with if = 2d, i.e. 

x,W = ix,(t)+ ^ (21) 

for some convenient vectors x^ ^ G ker (T'o (/?))• CH 

Let us consider first the well-posedness of the Dirichlet problem for C.E.L., 
i.e. (fT2|) and (fT3|) . We use some facts mentioned in 9, Section 3] which are 
consequences of the existence of the Smith form for regular QEP, see also [1, 

Let v — or v = n. Since admits exactly 2d distinct roots, then 

dimker(7'^(A)) = 1, VA G Q,, and the union of the various ker(7'i,(A)) spans 
C''. If jy — 71, we may decompose Xs(to) and Xs(ij) on the family {x^ q,}q^q^ 
and we obtain a linear system of 2c? equations w.r.t the 2d unknowns which are 
abscissas of Xs,a along the linear straight lines ker(7^„(a)), a G Q„. We proceed 
in a similar way when i> ^ and j G {1, . . . , n}. We get a linear system of 
2d equations by decomposing Xj(to) and JCj{tf) w.r.t. the 2d unknowns which 
are the abscissas of x^^, /3 G Qq. Hence, the Dirichlet problem amounts to 
solving n + 1 uncoupled square systems of size 2d (the very last one being 
useless due to the definition of x^). Each of the previous system is Cramer 
for almost all couple {tQ,tf). Indeed, the determinant of each system has the 
shape P{e^^° , e^*f )a6Q„ , v — or v ^ n, P being a polynomial with coefHcients 
depending on the coordinates of Xg^A and Xj^A- 

3.3. The case of D.E.L. 

Let us extend the Proposition 13.11 to the D.E.L. case. It should be empha- 
sized here that D.E.L. do not admit in general a unique solution. Nevertheless, 
given a solution, there exists one and only one pseudo-periodic solution which 
agrees with the first one on a grid Q^- As well as the study of autonomous 
dynamic differential systems leads to QEP, the study of autonomous difference 
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equations leads to transcendental eigenvalue problem associated to the following 
complicated matrix 

75.(e, A) :=-(Ji + 2(^.-1) J3) ^^e''^ - 

-2N < k < 2N 
-N < i < N 
\k + i\ < N 



^5 J2 J(7fc-7-fc)e'=^"-(-/2 + 2(i/-l)J4). (22) 

Let us introduce the following subsets 

{a e C/det(^^(£, A)) = o}, yv e M. 
Proposition 3.2. We assume that 

IQnl = IQol =47Vd, Q„nQo = 

detiVnie, 0)) ^ and det(^o(e, 0)) ^ 0. (23) 



Then there exists solutions x^ and x^ to ^14\ l and U5\) respectively of the shape 
(0) inside the interval [to + 2N£,tf — 2Ne]. 

Proof. The two last conditions imply that ^ Q„. If we set C = e^*^, we 
see that the quantity ■p„(e, X)C^^ is a polynomial w.r.t. ( of degree 4iV. So the 
equation C^^'^ det(P„(£, A)) = gives rise to a polynomial equation w.r.t. C, of 
degree ANd. 

Computation of the l.h.s. of ([T4| is performed by using ([Ij and [S, Lemma 6.1]. 
We find 

J2 ^n+aiXi{t)X-k{t){Ji+2in-l)J3)xs{t+ke) + iJ2+2{n-l)Ji)is{t) 

-2N < k < 2N 
-N < i < N 
\k + l\ < N 



^ 1 

+ y"\-k{t)-i-fk-l-k)J5is{t + ke)+nn^AJe + nJ7 = 0. (24) 

k=-N 



When t lies in the interval [to + 2Ne,tf ~ 2Ne], the various characteristic func- 
tions Xk{t) occuring in p4| are equal to 1. Next, we define for tj G [^o,*/] 
the grid Gtj,e = {tj + me,m e N} n [toj*/]- So, both restrictions of Xs(i) and 
Xj(i) to Qtj,e are vector- valued sequences satisfying linear constant matricial 
recurrences. The classical theory of those systems shows that, provided 

the characteristic equation admits a number of roots equal to the order of the 
recurrence, Xs(t) has the shape Xs(t) — Xs,o + ^\^'^*^s,\, Vt G Qtj,e, for some 
vectors x^ ,\ and x^.o- Here, the order of recurrence is equal to ANd and it is 
also equal to the number of roots of the characteristic equation which is | Qn | . 
So we may plug the previous formula into (j24p and we find 

^e^*^„(£,A)ic,,A+^„(e,0)i,,o =nn_elJ6 + nJ7. (25) 

A 
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Because the values of the function e , on the grid Gtj.e, are the numbers e ^^C™ 
with m = S N, all the functions e'*'' on this grid are linearly indepen- 

dent. Indeed, a linear relationship between these functions would give rise 
to a Vandermonde determinant w.r.t. to the associated distinct numbers C- 
Therefore, every non-constant function of t must vanish in (|25p . which means 
that the "phases" A occuring in Xs{t) are exactly the roots a of Q„. By as- 
sumption, 7^„(e, 0) is invertible and 7^„(e, a) is singular. Thus, we may choose 
Xs,o = n'P„(e, 0)~^(J7 -I- D-eUe) and Xs^\ € ker('P„(£, A)). Finally, we have 
determined the general solution Xg to (fT4|) on the grid Gtj ,s , namely 

i,(<) =x,,o+ ^"'^s^a- (26) 

Now, let us deal with ij{t). This function satisfies the following functional 
equation, which is similar to ([M)) 

X] ^n+am{t)x-k{t){Ji ~ '2.h)ij{t + ke) 

-2N < k < 2N 
-N < i < N 
\k + i\ < N 

^ 1 

-(J2 - 2J4)ij(<) - J5 ^x-fc W-(7fc - 7-fe)xj(t + ke) = 

2J3n_eais(i) + 2J4x,(t) - D-eUe + J7. (27) 

Let us construct a particular solution to (|27|) for t e By using the previous 

expression for Ks{t), the r.h.s. of (|27)) may be rewritten as 

□ _elJ6 - J7 - 2(n_eaiJ3 + J4)X.,0 - 2 X e"*( J3 + 0a J4)Xs,a 

where 6*0, = e~°'*D^^D^e"'^ — -p- X]fc j 7fc+ j 7j ^''"'^ • Now, if we substitute Xj{t) = 
Xj.o + X] e"*Xj.Q in ((27| . we note that the l.h.s. of (|27|) is equal to 

Po(e,0)%o+ e"*^o(£,a)xj-a. 

Because det(7^o(e: 0)) ^ 0, we may define 

%o = ^o(e, 0)-i(2(n_,aiJ3 + J4)xs,o - Q-sl + J7). 

Since Q„ H Qo = 0, the matrix VQ{e,a) is invertible for each a e Q„ and we 
may set 
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Similarly to the case of C.E.L., we readily prove that Scj^ — ;^Xs,o and x^ q = 
i-Xs^a, Va G Qn- At last, since |Qol = iNd, we conclude that the general 
solution to (P7)) on the grid Gtj ,e is given by 

x,W = lic,(t)+ ^ e^%^, (28) 

/3GQo 

for some convenient vectors x^ ^ in ker('Po(£, /?)). 

If we drop the requirement that t lies in Qtj,ei i-C- if we remove the con- 
dition ^-j^ e N, the functions t n- Xi(/:) and t ^ Xj(t) may be extended by 
the preceding formulas to pseudo-periodic functions t i-^ Xs(0 and t i-^ ij{t) 
respectively. Since the equations of motion are autonomous (independent w.r.t. 
t), these n+1 functions are solutions to ([M)) and ([Tsl) respetively. Therefore, 
these functions are of the shape ([3]) with K ~ ANd and K = 8Nd respectively 
and the proof is complete. □ 

Remark 3.1. Solving D.E.L. with Dirichlet conditions leads to n -f 1 uncoupled 
linear systems, one of size X^aeg dimker(7^^(£, a)) and the n others of size 

J2i3eQ„ dimker(7^o(e, Z?))- If those systems are Cramer, then the pseudo-periodic 
solution to D.E.L. exists and is unique. 



4. Convergence issues 

Let us fix jy, 7_7v, . . . , jpf, Ji, . . . , J5. Motivated by studying the convergence 
of the solutions to D.E.L. to the respective solutions to C.E.L., it is natural at 
first sight to ask if the matrix- valued function Vi,{e,X) tends to Vu{X) locally 
uniformly w.r.t. AsCase tends to 0. Next, we recall the Hausdorff metric 

dniFi, F2) '■= max < max minja; — y], max min — ?/| >, 

defined for all nonempty finite subsets -F'i,i^2 C C. Thus, we naturally investi- 
gate the convergence, in this sense, of 

= (det(^,(e,.)))-i{0} to = {Aet{V.{.)))-^ {Q} 

as £ tends to 0. In order to prove this result, we shall need the following Theorem 
of Cucker and Corbalan l]. 

Theorem A.l. Let P{X) = aoX™+aiX™-i + . . .+«,„ G C[X]\{0}. Let^i,...,^r 
be its roots in C, with multiplicities Hi, . . . , fir respectively, and let Bi, . . . ,Br be 
disjoint disks centered at ^1 , . . . , ^j. with radii Eq and contained in the open disk 
centered at with radius I/sq- Then, there is a S £ R^, such that, if \bj — aj\ < S 
for every < j < m, then the polynomial Q{X) = boX"^ + biX™~^ + . . . + bm 
has fj,j roots (counted with multiplicity) in each Bj and Aeg{Q) — deg(P) roots 
with absolute value greater than I/eq- 
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It extends older results of Weber and Ostrowski to the case of perturbation of 
polynomials of distinct degrees. Hence, we must exclude the 4Nd— 2d divergent 
roots, as e tends to 0, from the set Qi, to prove the second result of convergence 
mentioned above. 

Theorem 4.2. We keep the assumptions il9\) and \23l\) of Provositions \3. 1\ and 
We assume that defined by (QP is such that 7-jv7jv 7^ and 

□jj , e [to + 2Ne,tf ~ 2Ne]. (29) 

Let v Cz R. and K any compact neighbourhood of . Then, when e tends to 0, 
'Py(e,A) tends to 'P^{\) locally uniformly in (C and Qi/ H K tends to Qjy in the 
Hausdorff sense. 

Proof. The assumptions (j29p are equivalent to the algebraic equations 

N I ^ 

^ 7fc = and - ^ k{jk - 1-k) ^ 1 

k=-N k=~N 

since the characteristic functions are equal to 1 in [to + 2Ne,tf — 2Ne]. In The- 
orem 6.1 of fs], we have proved that these conditions are themselves equivalent 
to one or the other statements 



for all X G C'^ {[to , t f]) , lim 0^x{t) = x(i) locally uniformly in ]to,t 



fl 



• for all x e C'^ {[to , t f]) , lim 0^^x{t) — — x(i) locally uniformly in ]to,tf[. 

The mode of convergence means that for all d > 0, □±ex(t) tends uniformly to 
±x(t) in [to + <5, i/ — S] when e tends to 0. This convergence can not be improved 
since the functions 1 1-^ D^l and 1 1-^ D^t are equal to and 1 respectively only 
in the interval [to + 2Ne,tf — 2Ne]. By composition of these properties we 
obtain 

- = - E ^ik+aixAt)x-k{t)e'''' -f^x^ (30) 

-2N < k < 2N 

-N < ^ < N 
[k + e[ < N 



and 



^(□ee^*-n_,e^*)= f2xMt)l{lk~-l-k)e''^^^2\. (31) 

k=-N ^ 



We see easily that the functions e~^*'D^^\J^e^^ and e~'^*(nee^* — □_ee'^*) defined 
at e = by the respective values — and 2A are continuous w.r.t. e. The 
quantities in both sides in each equation are obviously the coefficients of ( Ji + 
2(i/-l)J3) and J5 in ([TH) and (HU when i lies in the interval [to+2Ne,tf~2Ne]. 
Hence, the mapping A Vp{e, A) tends to A i— )• "Pi/ (A) uniformly on any compact 
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subset of the preceding product, as e tends to 0. 
Let us deal with Q^. We compute first 

V,[e, A) - V,{\) = ( Ji + 2{u - l).h){\^ + e-^*n_,ae^*) + 

J5(-2A + e-^'U,e^' - e-^'U^,e^'), (32) 

for all e ^ 0, t G [^o,^/] and A G C. The l.h.s. of ([5^ is independent of 
t and the r.h.s. is constant w.r.t. t inside [ig + 2Ne,tf — 2Ne], as we see 
from ([50)1 and ([5T|) . Expanding in Taylor series the exponentials e*^'^^ w.r.t. 
e, we find a matrix-valued convergent Taylor series w.r.t. e for ^'^(e, A). The 
coefficient of e™ in V^{e, A) is a polynomial matrix w.r.t. A, independent of t 
inside [to + 2Ne,tf — 2Ne]. Now, the determinant of such a convergent Taylor 
series is itself a convergent Taylor series. 

At this point we have established that ^^wd (Jet(7^y (e, A)) is a polynomial of 
degree 4Nd w.r.t. = e^'^ and admits a Taylor expansion w.r.t. e starting at 
C2^'^det(7'.(A)). 

We choose Eq so that Q,y C K C K C B{0, I/sq) and small enough to separate 
the elements of Q^. Let S as in Theorem 14.11 We choose next e so that, if 
1 < TO < iNd, the coefBcicnt of C™ in det{Vu) - det('P^) is less than S. Now, 
we may formulate the conclusion of Theorem 14.11 as the following inclusion 

Q,-(det(7^,(£,.)))-H0}c (C\S(0,l/£o))U [ |J 6(A,eo)). 

As a consequence, intersecting both sides with K we get dniQu ^ K,Qu) < eq 
for all £ small enough. This ends the proof. □ 

Remark 4.1. The convergence of x(t) to x(i) as e tends to implies more 
complicated issues. Indeed, not only the phases Q„ and Qo have to tend to 
Qn and Qo respectively but the amplitudes Xs,q, Xj^q, and x^-^, where a € Q„ 

and (3 G Qo, have to tend also to the respective amplitudes Xs,a, Xj,a, and x^^^, 
where a € Q„ and /3 € Qq. We refer to [Sj for an examination of the difficulties 
in the case n = 1. 



5. Periodicity and choreographies 

We focus in this section on periodic and choreographic solutions. Let us 
define a choreography of n particles (xi, . . . , x„) in R'* as a T-periodic solution 
to the equations of motion in which the trajectories differ one to the other by 

kT 

some delay of the shape — , k G {1, . . . ,n}. In other words, a choreographic 

n 

solution is a mapping u : R/TZ — > R'^ such that u{t + T) — u(t) and such 
that the family {xj{t)}j, defined by Xj(t) ~ u(t + jT/n), satisfies for all j 

u{t+^) - F,(u(t+^),...,u(t + r)), 
-□_.au(t+f) - F,(u(t + ^),...,u(i + r)), 
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i.e. the respective equations of motion C.E.L. and D.E.L. presented in ([5]). 

Theorem 5.1. 1. Under the assumptions of Provosition [K7[ all the solutions 
to C.E.L. are periodic if and only if 

Qn U Qo C iR, VA', A" e Q„ U Qo, A'/A" e Q, (33) 

2. Under the assumptions of Provosition WM all the pseudo-periodic solutions 
to D.E.L. are periodic if and only if 

Qn U Qo C VA', A" e Q„ U Qo, A'/A" e Q. (34) 

3. //det(7'„(0)) ^ 0, det(-p„(e,0)) ^ and 

IQol = 2d, Qo C VA',A" e Qo, A'/A" e Q, (35) 
IQol = 4iVd, Qo C iR\ VA',A" e Qo, A'/A" e Q, (36) 

f/ien there exists choreographic solutions Xj(t) and Xj(i) to C.E.L. and 
D.E.L.. 

Proof. 1. We first notice that if {ue}e=i,....,K, is a family of nonzero vectors in 
C^, then the various functions 1 1— >■ e'''**U£ are linearly independent iff the 
Xe are pairwise distinct. It relies on the nonsingularity of the Vandermonde 
matrix V(Ai, . . . , Xk). As a consequence, the function u{t) = J2f=i e'^**Uf 
is periodic iff for some T > we have Vi, XiT e 2i7rZ and this is equivalent 
to the requirement V€, iXi G M* and Vj, fc, Xj/Xk & Q*. Therefore, the 
period T of u{t) is inf{T > 0, G Z,V^}. Taking in account that 
the vectors Xs^„ and x^-^, occuring in the proof of Proposition 13.11 may 
be chosen arbitrarily in the respective appropriate null spaces ker('P„(a)) 
and ker(7'o(/5)), the previous properties of periodicity apply to the set of 
solutions to C.E.L. and give formula ((33)) . 

2. Pseudo-periodic solutions to D.E.L. are of the shape ([3]) by using Propo- 
sition [32] and the previous arguments apply. 

3. We first observe that for each choreographic solution of the shape ([3]), 
Xs(t) is necessarily constant. Indeed, 

J = l V " / £=1 j = l 

By periodicity, we have XiT e 2i7rZ for all £ so that }J}Li " j = 0. 

Having this fact in mind, we may solve (|13p . Our assumptions imply 
that the solution Xj(t) to (IT^ may be written as ([3]) with J-iT = 2(i, 
since the underlying Quadratic Eigenvalue Problem satisfies |Qo| — 2d 
and det(7'„(0)) ^ (see Proposition [3l]) . Plugging Xj(t) = u(i + ^) 
into (|13p and using the linear independence of the summands ([3]) , we see 
that is satisfied if and only if Uo = 7'„(0)~^ J7 and for all A € Qo, 
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u\ G ker('Po(^))- If we choose the vectors Xj{to) and Xj(to) or Xj(to) and 
'x.jitf) for all j according to the preceding explicit form for Xj(t), we have 
justified the existence of choreographic solutions to C.E.L. with Xs = est. 
Let us deal now with D.E.L.. As seen in the proof of Proposition 13.21 
each solution Xj{t) to has the shape ([3]) with K = ANd due to our 
assumptions on Qq and 'P„(e, 0). The remainder of the proof is entirely 
similar to C.E.L. . First, ([TS]) is satisfied if and only if Uq = Vn{£, 0)~^{J7 — 
□ -eJg) and for all A G Qq, u\ £ keT{Vo{e, X)). Second, convenient choice 
of initial or boundary conditions guarantee the existence of choreographic 
solutions to D.E.L. with x^ — est. 

□ 

Remark 5.1. We may convert the existence of choreographic solutions into a 
linear algebra problem. Indeed, we add to the systems described at the end of 
the Sections 3.1 and 3.2 the following equations 

Xs,a = 0, Xj,Q = e"(^~i)n xi,a and x^-^^ = e^(^~i)ixi,,3, 

for all a £ Qn,P & Qo and j = 1, . . . ,n. Due to (^0)) and (PTt we see that, 
provided Dirichlet problem is well-posed, we find a choreographic solution. 

Remark 5.2. In the litterature (see for instance j§. Section 3.10]), the problem 
of the existence of choreographic solutions arises when one studies gyroscopic 
systems. The algebraic conditions P7 and P8 in Table 1.1] amount to require 
that Ji = Vi > 0, J2 = V2 > 0, J3 and J4 symmetric and small enough 
compared to Ji and J2 respectively. 



6. Numerical experiments on choreographies 

Experimental and working algorithms performed in this last section are 
implemented in Maple and Matlab. We deal with real symmetric matrices 
Ji, . . . , J4, zero vectors Jg, J7, small dimension systems {d — 2, 3) since it dis- 
plays already the main features, and arbitrary number of particles. Furthermore, 
existence of periodic or choreographic solutions requires that J5 = 0. Let us give 
some details on the choice of the matrices Ji. Given Ji, J2, J3, we set, if d = 2, 

J4 — i J2 -I- i( Ji — 2J3) ( ■^.^ ■^.^ I . Identifying the coefficients of the polynomial 
\J3 .HJ 

det(7'o(A)) with those of (A^ -I- /?i)(A2 + /3|) and requiring that J4 = V4, we 
get three equations on ji, J2, J3, the coefficient j'4 standing free. Thus, we may 
choose Ji -|-2(n— 1) J3 and Ji — 2J3 definite positive, J2 — 2J4 and J2 -l-2(n— 1) J4 
definite negative. 

We present in Figures [T] and [5] the graphs of two solutions to gyroscopic 
C.E.L., sharing the same matrices Ji, J2, J3. On the left, a typical periodic curve 
obtained by considering [3i = Ai and (32 = lOi and on the right, a non-periodic 
curve. Incommensurability between Ai and 7i^/2 explains the non-choreographic 
behaviour of the curve, as mentioned in property p3|) . 
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Figure 1: C.E.L., /3i = 4i and 132 = Wi 



Figure 2: C.E.L., /3i = 4i and = YiV^ 



Let US deal now with D.E.L.. For sake of clarity, we shall denote by Xj{t) 
and yj,Af(i), Vj e {1, . . . ,n}, the unique solution to C.E.L ([T3| and the unique 
pseudo-periodic extension to [Lq, tf] of the unique solution to D.E.L. (j27p on the 
grid tJto,e with e = ^^^^g^ respectively. 

First, we give some hints to solve (|27l) . When t e [to + SiVe, t/ — 2iVe], we may 
compute yj{to + 2Ne) as a function of yj(io + ke) with fc varying from —2N to 
2 — 1 . If t ^ [to + 2Ne ,tf — 2Ne] , some of the characteristic functions occuring 
in (|27|) vanish and solving (f27| must be slightly modified, see more details in 

[§| . We consider the matrices Ji — ^ , = ^ , J3 = 

and (/3i,/32) = (2i,5i). In a first experiment, we use an operator De^fc such 
that = 1 and (7-1,70,71) = (— ^ + ik,—2ik, i + ik) where k £ R. Because 
Xj+i(t) — Xj{t + iT) and yj+i,A/(t) = yj,M{t + -^T), all the particles have the 
same trajectory, either in both cases C.E.L. and D.E.L. . Figure [3] depicts the 
curves of yj-,35, yj,42, yj,75 and x^ . 

Since our algorithms are suitable for each dimension, we provide also an 
example of quadratic choreography with d = 3 in FigureHl We choose Ji, . . . , J4 
such that det(Po(£, A)) = (A^ + l^A^ + 4)(A2 + 9) and the same operator D^^k 
than previously. 

Figures [3] and |4] illustrate roughly the phenomenon of convergence, as M 
increases, of the pseudo-periodic solution yj^M{t) to D.E.L. to the solution Xj{t) 
to C.E.L., for all j, for all operator De.fc. In order that the approximation 
becomes meaningful, the number AI — ^[Qi^ ,,) — 1 must satisfy 

M2 » ^^-^{tf - hfp{{Ji ~ 2J3)"i(J2 - 2J4)) 
l7-i7i| 

as one sees from equation (P7| (p denotes here the spectral radius). Let us 
remark that this lower bound is independent on d and still holds for general 
operators 

In the next experiment, we still work with the previous lagrangian, with 
= 1 and with an operator such that 70 = — (7-1+71 )• We choose M = 100 
in order to avoid some erratic behaviour observed, for example, in the first plot 
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Figure 3: Solution of D.E.L for M = 35, 42, 75 and solution to C.E.L. in 




Figure 4: Solution of D.E.L. for M = 30 and solution to C.E.L. in 



(M = 35) of FigureEl We provide the plot of - log(niin(||xj (t) - yj(i)||2, 3Af)) 
as a function of (7_i,7i) G in Figure [5j Two peaks occur at (7_i,7i) = 
±(1/2,1/2) and reveal a good approximation of Xj (t) by (t) . The two previous 
pairs are better understood if we have a look to the error 'x.j{t) — yj{t) with 
complex operators iDg^fc. In Figure [51 we give the plot of the 2-norm of the 
error with (7_i, 70, 71) = (— 5 + ik, —2ik, ^ + ik) for several values of M. The 
operator De^o — seems to be in any case the better choice. 

Let us conclude this paper with an additional remark on the convergence of 
solutions, completing Remark l4.ll Recall first that D.E.L. converges to C.E.L 
iff if of the shape ([1]) and checks D^l = and Dei = 1, provided t G 
[to + 2Ne,tf — 2Ne], see 0, Definition 6.1. and Theorem 6.3]. In that case, 
the condition (7-1,70,71) = (—5 + ik, —2ik, ^ + ik) is linked to the inclusion 
Qo C jM* as mentioned in [1, Proposition 5.1]) for the special case d = 1. Based 
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Figure 5: - log(mm(||x^ (t) - y, {t)||2 , 3M)) ^.^^^^ ^^^^^ ^^^^^ ^.^^ ^ 

as function of 7-1 and 71 (M = 100) ^^^^^^j ^^^^^^ ^ 



on the preceding experiments, we conjecture that under mild condition of non- 
resonance of the lagrangian, the solution to D.E.L. converges to the solution to 
C.E.L., as e tends to 0. 
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